Relativistic coupled-cluster single-double calculations of positron-atom bound states 
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Relativistic coupled-cluster single-double approximation is used to calculate positron-atom bound 
states. The method is tested on closed-shell atoms such as Be, Mg, Ca, Zn, Cd, and Hg where a 
number of accurate calculations is available. It is then used to calculate positron binding energies 
for a range of open-shell transition metal atoms from Sc to Cu, from Y to Pd, and from Lu to Pt. 
These systems possess Feshbach resonances, which can be used to search for positron-atom binding 
experimentally through resonant annihilation or scattering. 

PACS numbers: 36.10.-k, 34.80.Uv, 34.80. Lx, 78.70.Bj 



I. INTRODUCTION 

In this paper we apply the relativistic coupled-cluster 
single-double (SD) approximation to calculate positron 
binding energies for a number of atoms, including open- 
shell transition metal atoms. Positron-atom bound states 
are characterized by strong electron-positron correlation 
effects. These effects for the positron bound states are 
even stronger than for their electron counterparts, i.e., 
the atomic negative ions. This is due to the electron- 
positron attraction and, in particular the large role of 
virtual positronium (Ps) formation I1H3|- This makes 
calculations of positron-atom bound states a challenging 
theoretical problem. Using the SD method allows us to 
include all the main correlation effects for atoms with 
several electrons in the open valence shell, which has not 
been done before. 

The existence of positron-atom bound states was pre- 
dicted by many-body theory calculations Q and veri- 
fied variationally [g, |6[ more than a decade ago. Since 
then positron binding energies have been calculated for 
many ground-state and excited atoms: He 2 3 S, Li, Be, 
Be 2 3 R Na, Mg, Ca, Cu, Zn, Sr, Ag, Cd and Hg (see, 
Refs. [7Hl2| and references therein). There are strong 
indications that many more atoms, possibly over a quar- 
ter of the whole periodic table, should be able to bind 
positrons 0, [j| [lj, [l4| • Such conclusions are based on 
the analysis of the atomic ionization potentials and dipole 
polarizabilities, which are relevant for positron binding, 
and on the calculations of positron binding to a "model 
alkali atom" (TBj (see also review on positron compounds 
,16j). In spite of this wealth of predictions, experimental 
verification of positron binding to neutral atoms is still 
lacking. 

Recently we proposed (Ref. [l3j]) that positron bind- 
ing to many open-shell atoms could be studied experi- 
mentally by measuring resonant positron-atom annihila- 
tion. Such resonant annihilation should be similar to that 
observed for positrons in many polyatomic molecular 
species [I?) ■ In this process the incident positron is cap- 
tured into the bound state with the target, with the ex- 
cess energy being transferred (in the case of molecules) to 



vibrations. Since the vibrational motion of the molecules 
is quantized, these transitions can only take place at spe- 
cific positron energies, which means that they have a 
resonant character. These energies correspond to vibra- 
tional Feshbach resonances of the positron-molecule com- 
plex The majority of the resonances observed 
are associated with individual vibrational modes of the 
molecule. The negative energy of the positron bound 
state eo is then related the downshift of the resonance 
energy, 



eo, 



(1) 



with respect to the energy io v of the vibrational excita- 
tion [2(| Hence, by observing the resonances, the 
positron binding energy Sb — |eo| can be found. In this 
way binding energies for over sixty polyatomic species 
have been determined [22|-|25| by measuring positron an- 
nihilation using a high-resolution, tunable, trap-based 
positron beam [26| . In our previous paper [l3| we sug- 
gested that one can search for a similar effect in atoms. In 
this case the resonances will be associated with low-lying 
electronic excitations. These can be found in open-shell 
atoms, where they often have the same configuration as 
the ground state. If the positron can bind to such an 
atom in the ground state, then it is likely to bind to its 
excited state as well. One can then consider the following 
process, 
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Here the positron first loses energy by exciting the atom 
and becoming trapped in the bound state with the ex- 
cited atom. It then annihilates with one of the electrons, 
and the resulting gamma quanta can be detected. The 
first step of the process (O is obviously reversible. Hence, 
to estimate the efficiency of resonant annihilation one 
needs to evaluate the rates of both positron annihilation 
and autodetachment [HI]. 

To prove the feasibility of such process one also needs 
to estimate the positron binding energies to open-shell 
atoms. In Ref. |13| we used a simple approach in which 
the second-order correlation potential of the positron- 
atom interaction was scaled up to account for the effect 
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of higher-order electron-positron correlations (i.e., virtual 
Ps formation). It was assumed that the scaling factor is 
the same for all atoms, and its value was chosen by fitting 
the results for those atoms where accurate positron bind- 
ing calculations are available. Open shells were treated 
by using fractional occupation numbers in the standard 
expressions for a closed-shell system. This approximation 
looks reasonable in the positron-atom problem, because 
the positron is not affected by the Pauli principle. The 
positron-atom interaction has, thus, no direct sensitiv- 
ity to the valence shell being open or closed (unlike the 
electron-atom interaction in the negative ion problem). 

In the present paper we use a more sophisticated fully 
ab initio approach in which the strong electron-positron 
correlations are included explicitly through the use of an 
all-order technique. 

Accurate treatment of the strong electron-positron cor- 
relations is the main challenge in the calculations of 
positron-atom interaction. It calls for the use of non- 
perturbative approaches. For systems with few active 
electrons accurate results were obtained using stochas- 
tic variation method [fl Q . The most obvious choice for 
a generic many-electron atom is the configuration inter- 
action (CI) technique. It was successfully employed in a 
number of previous calculations for atoms with one or two 
electrons in the valence s shell (see, e.g., @, H, H3, HI]). 
However, it becomes too complicated for more than two 
valence electrons and, in particular, for atoms with open 
d-shells, which we want to consider in present work. An- 
other suitable all-order technique is the couple-cluster ap- 
proach. In this approach the interparticle interaction is 
included to all orders via an iterative procedure. The 
corresponding subset of terms includes the so-called lad- 
der diagrams [29] . This class of diagrams is very im- 
portant in the positron-atom problem, as it describes the 
effect of virtual Ps formation. Summation of the electron- 
positron ladder-diagram series was performed earlier by 
solving a linear matrix equation for the electron-positron 
vertex function for hydrogen [3(|, noble-gas atoms [3lj . 
and halogen negative ions (32| . 

The coupled-cluster approach in its single-double (SD) 
approximation has been widely used for atoms and ions 
with one external electron above closed shells (see, e.g., 
[33M37I ] ^ . It is relatively easy to modify the corresponding 
equations for the case of a positron interacting with a 
closed-shell atom. We will do this in the next section. 

To test the method we first apply it to positron bind- 
ing to closed-shell atoms: Be, Mg, Ca, Zn, Cd, and Hg 
for which a number of accurate calculations is available. 
We then apply the same approach to open-shell atoms, 
treating them in a simplified manner similar to that of 
Ref. [131 ] . Fractional occupation numbers are used to 
rescale the terms containing contributions from the open 
shells. The main advance of the present method com- 
pared with Ref. [IH is the all-order treatment of the 
positron-electron interaction in the correlation potential. 
No further fitting is used or needed. Remarkably, the re- 
sults of the all-order calculations turn out to be close to 



the estimates obtained in our previous work. These cal- 
culations lend further support to the proposal to search 
for positron-bound states with open-shell atoms through 
resonant annihilation or scattering. 



II. THEORY 

A. Electron and positron basis sets 

The use of the coupled cluster technique requires con- 
struction of a single-particle basis. For a positron inter- 
acting with an atom one must have two separate basis 
sets - one for electron states and one for positron states. 
We use a standard B-spline technique in both cases [38[ . 
The electron basis states are constructed by diagonal- 
izing the matrix of the relativistic Hartree-Fock (RHF) 
Hamiltonian in the B-spline basis. The positron basis 
states are constructed using the same set of B-splines 
and the RHF Hamiltonian in which the sign of the direct 
potential is changed and the exchange potential is omit- 
ted. Below we will use the following notation for the basis 
states: indices o, b, c refer to electron states in the core, 
indices m,n, k,l refer to electron states above the core, 
indices v, r, w refer to positron states, and indices i,j re- 
fer to any states. Numerical results reported in Sees. IIIII 
and IIVI were obtained using the basis sets built from 30 
B-splines of order 7 spanning the radial coordinate from 
the origin to R = 30 a.u., with angular momenta between 
and 10. 



B. Singles-doubles equations 

The wave function of an atom with a positron in state v 
can be written in the SD approximation as an expansion 
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where is the zeroth-order wave function of the 

frozen-core atom in the relativistic Hartree-Fock approx- 
imation with the positron in state v. It can be written 
as 



= at 
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(4) 



where |0c) is the RHF wave function of the atomic core. 

The expansion coefficients p na and p mn ab in Eq. © 
represent single- and double-electron excitations from the 
core. The coefficients p rv represent excitations of the 
positron, and the coefficients p rnW a represent simultane- 
ous excitatons of the positron and one of the electrons. 

The SD equations for the core excitation coefficients do 
not depend on the external particle and are the same in 
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the electron and positron cases. They are written as a set 
of equations for the single-excitation coefficients p m „ and 
double-excitation coefficients p m nab (see, e.g., Ref. }33|). 
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and 



(6) 



(^a H~ ^n) Pmnab — Qmnab 
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cd kl 
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In these equation <7 m „fc/ = g mnfc / - 5 mmtt and /5 m „/d = 
Pmnki — Pmnik, and the coefficients g are the Coulomb 
integrals, 

9mnki = ( /^(ri)^(r 2 ), — r^fc(ri)V'i( I '2)dr 1 dr2, 

77 |r x — r 2 | 

(7) 

involving electron states ipk, ^Pu etc. 

The core SD equations and © are solved iteratively 
until convergence is achieved. This is controlled by the 
correlation correction to the energy of the core, 
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After solving the SD equations for the core one can 
start iterating the SD equations for the external parti- 
cle. The SD equations for the positron can be obtained 
by substituting the state \^ v ) from Eq. ([3]) into the rel- 
ativistic many-body Schrodinger equation, 



H\V V ) =e |*„>. 



(9) 



Projecting the Schrodinger equation onto aJjOc) gives 
the equation for p wv , 

(^0 ^w)Pwv — ^ QwbvjnPmb ~l~ ^ QwbrmPrmvb: (1 0) 



Projecting Eq. (|9]) onto a^a^aa^c) gives the equation 
for the double-excitation coefficient p wri va, 
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FIG. 1: Diagrammatic form of the SD equation (|1U|) for the 
positron single-excitation coefficient. Lines with arrows to the 
right are positron and excited electron states, and those with 
arrows to the left are core electron states ("holes"). The verti- 
cal solid line is the double-excitation coefficient, and the solid 
line terminated by a circle is the single-excitation coefficient. 
Dashed lines are the Coulomb interactions. 



Here q wnva is the Coulomb integral ([7]) involving positron 
states. 

Equations (| 10[) and (TTTj) are presented graphically in 
Figs. [T]and[2J When solving these equations, the correc- 
tion to the energy of the positron state v, 



fi^v ^ ^ QvbvmPinb ~t~ ^ ^ QvbrmPrmvbi 

brar 



(12) 



nib 



is used to control convergence. 



C. Positron bound state 

Apart from the absence of the exchange terms involv- 
ing the positron, there is another important difference 
between the SD equation (fnj) for the positron and those 
for external electrons in atoms or ions. In the latter cases 
the SD equation for the double-excitation coefficient also 
contains a term ^ r q wnra Prv The corresponding dia- 
gram is "reducible", i.e., it contains two parts connected 
by the single line of the excited valence particle. Includ- 
ing it ensures that the correction to the energy of the 
valence state v is determined in all orders. 

When using the SD method for the positron, its eigen- 
state and energy eigenvalue eo are obtained by matrix 
diagonalization (see below). The need for this procedure 
arises because of the absence of a good zeroth-order ap- 
proximation for the wave function of the bound positron. 
In the RHF approximation the positron-atom interac- 
tion is repulsive, and all of the single-particle positron 
basis states lie in the continuum. On the other hand, 
the B-spline basis is effectively complete for the positron 
bound state to be obtained by matrix diagonalization. 
To do this, the SD equations (TIT))) and pip should be 
iterated for every positron basis state used to construct 
effective single-particle Hamiltonian matrix. A similar 
situation occurs for negative ions [3 3 . |4fj and when com- 
bining the SD method with the configuration interaction 
technique to obtain many-electron wave functions [4l| . 
In all of these cases the term J2 r IwnraPrv should be re- 
moved from the SD equations, as it is taken into account 
via matrix diagonalization. 
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FIG. 2: Diagrammatic form of the SD equation for the positron double-excitation coefficient. 



The wavefunction of the bound positron is found as an 
expansion over the positron basis states, 



D. Third-order corrections 



(13) 



Equations (|10l) and (ITT|) are iterated for every basis state 
v in the expansion (1131) . The energy parameter eo m these 
equations is the (unknown) energy of the positron state 
(|13jl . The energy eo and the expansion coefficients c„ are 
found by solving the eigenvalue problem 



El = e X, 



(14) 



where X is the vector of expansion coefficients c v , eo is 
the lowest eigenvalue (which must be negative), and the 
elements of the effective Hamiltonian matrix £ are given 
by 
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(15) 



The first term on the righ-hand side of Eq. (fT5]) repre- 
sents the positron energies in the static RHF approxima- 
tion. The second and third terms describe the effect of 
electron-positron correlations. These two terms are the 
positron-atom self-energy, which is given by right-hand 
side of the diagrammatic equation in Fig. [T] 

Since the SD equations (fTU|) and ([TT|) depend on the 
energy eo which is found later from Eq. fjl4[) . we start 
with an initial guess for eo. The calculations are then 
performed iteratively, solving the SD equations (fTOj) and 
(|11|) and diagonalizing the matrix f|15|) several times until 
eo has converged. In practice this takes about 5 global 
iterations. 



The SD equations account for the second-order contri- 
bution and selected classes of higher-order correlation di- 
agrams in all orders. In particular, the electron-positron 
ladder diagram series, which describes virtual Ps forma- 
tion, is included in full. However, beginning with the 
third order, the SD approximation misses certain terms. 
It is well known that the missing third-order terms can 
give sizeable corrections to the energy in atomic systems 
(see, e.g., Ref. 35]). Including these terms can lead to 
significant improvements in the accuracy of the results, 
and we include these contributions for the positron bound 
states with atoms. 



The third-order terms not included in the SD approx- 
imation are listed in Ref. [35| for the case of electron 
interacting with a closed-shell ion. Similar terms for the 
positron case can be obtained by changing the sign of 
Coulomb integrals involving the positron and removing 
exchange. All resulting third-order diagrams are shown 
in Fig. |31 The corresponding perturbation-theory correc- 
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FIG. 3: Third-order positron-atom interaction diagrams not included into the SD equations. Notation of Ref. [43 ] is used. 



tions to the energy of the positron state p are 
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Here = e, + e_y, eo is the positron energy found by 
solving the eigenvalue problem (]14p . and the positron 
state p is given by Eq. ([131) . 



III. CLOSED-SHELL ATOMS 

For a positron interacting with a closed-shell atom all 
sums over the projections of angular momenta of elec- 
tron and positron states in the SD equations ([5]), ([6]), 
(ITU! , and (fTTj) and third-order corrections (lll)l) - ([2"0"[) . are 
done analytically, together with the angular reduction of 
the Coulomb integrals. Such systems are easiest from 



the computational point of view. There are also a num- 
ber of other calculations of positron binding of different 
level of sophistication for several such atoms. We use 
these systems to test our approach. The results of calcu- 
lations of positron s-wave bound states with Be, Mg, Ca, 
Zn, Cd, and Hg are shown in Table Q] These six atoms 
share the following features. First, they have closed elec- 
tronic shells. Secondly, the dipole polarizabilities ay of 
these atoms are sufficiently large to provide strong attrac- 
tion and ensure positron binding. Finally, the ionization 
potentials I of these atoms (with the exception of cal- 
cium) are greater than the Ps binding energy (6.8 eV). 
Hence, the Ps-formation channel is closed at low positron 
energies, and the positron-atom bound states are stable 
against dissociation into positive ion + Ps. 

The SD column of Tabic fl] shows the binding energies 
(BE) | eo | (in meV) found by solving the eigenvalue prob- 
lem (ll~4"l) . The E3 column gives the third-order correc- 
tions (fT^[) ([2"U[) . The next column is the total obtained 
by adding the third-order correction to the BE from the 
SD calculations. We also examined the effect of finite 
box size on the binding energy [7j, and found it negli- 
gible. The rest of the Table shows the results of earlier 
calculations. The most accurate among them are proba- 
bly the calculations by the Mitroy group [27l. [28 ll43l - l45j . 
Our results are generally slightly larger, but within 100 
meV of the best earlier predictions. Adding third-order 
corrections increases the BE in all the systems, and in 
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TABLE I: Positron binding energies (BE, in meV) for closed-shell atoms obtained using the SD equations (SD) with third-order 
correction (E3), and by other methods, with best available predictions in bold. 



Atom 


I 


Old 


BE, this work 


BE from other calculations 




(eV) 


(a.u.) 


SD 


E3 


Total 




4 Be 


9.32 


38 


157 


26 


184 


45.90 [43], 85.63 [44], 83.8 [45], 3.0 [46], 33 [47] 


12 Mg 


7.64 


72 


475 


37 


512 


870 [4J, 425 [44], 439 [27], 985 [48], 15 [46], 457 [47], 125 [49] 


20 Ca 


6.11 


154 


1114 


50 


1164 


1139 6 [2JJ 


30 Zn 


9.39 


42 


143 


40 


183 


230 [4J, 39 [50], 103 [28], 0.01 [51], 53 [49J 


48 Cd 


8.99 


49 


204 


70 


274 


350 [4], 166 [28], 1.5 [51] 


80 Hg 


10.43 


38 


59 


53 


112 


45 [4] 


"Static dipole 


polarizabilitics from Ref. 


& 









c Positron-atom BE obtained from the calculated BE of 449 meV 
with respect to the Ca + +Ps threshold. 



most cases increases the difference between the present 
and best earlier results. 

Note the scatter of the results between different groups 
and even between different calculations by the same 
group. This reflects the fact that the positron- atom bind- 
ing energy is very sensitive to the correlations. In the 
static field RHF approximation the positron is repelled 
from the neutral atom. The binding is solely due to cor- 
relations. Most of the correlation energy is required to 
get the system across the threshold for binding, while the 
final binding energy is a result of a relatively small "sur- 
plus" of the correlation energy. Note that according to 
our calculations, positrons do not bind to either Pd or Pt 
in the atomic ground states. However, both atoms bind 
in excited states. 



IV. OPEN-SHELL ATOMS 

In this section we consider open-shell atoms which were 
suggested in our previous work [l3| as good candidates for 
experimental detection of positron-atom bound states via 
resonant annihilation or scattering. We treat open-shell 
systems in an approximate way, by introducing fractional 
occupation numbers. For example, the ground-state elec- 
tron configuration of neutral iron is 3d 6 4s 2 above the 
Ar-like core. We treat it as a closed-shell system but 
reduce the contribution of the 3d subshell to the cor- 
relation potential by the factor 0.6. Both members of 
the fine-structure multiplet, 3g? 3 / 2 and 3c?5/2, are treated 
identically and the corresponding terms are rescaled by 
the same factor 0.6. 

Note, however, that the SD equations (fTTj) are left un- 
changed. There is strong cancellation between different 
terms in these equations and rescaling of the core contri- 
bution is an insufficiently accurate procedure. It breaks 
the delicate balance between different terms, leading to 
unreliable results. The rescaling is done when construct- 
ing the correlation potential matrix a vw , by reducing the 
terms corresponding to the open subshells when summing 
over the hole states b in Eq. (IT5|) . The fact that this 
procedure works well is supported by good agreement 



between the present calculations of positron binding to 
copper with our previous calculations (see Table [TTJ) . 

This approach does not distinguish between different 
states of the same configuration. This means that if the 
positron is found to be bound to an atom in a given 
electronic configuration, it is predicted to be bound to all 
states of this configuration. This should be expected at 
least for the states where the excitation energy is smaller 
than the binding energy. However, this may still be true 
for higher-lying states. The ability of an atom in an 
excited state to bind the positron is what is needed for 
the resonant process proposed in Ref. [l3j |. As it was 
argued in Ref. [13| and as we will see below, there are 
many such atoms. 

Table [XT] shows the results of our calculations of 
positron binding energies for a number of atoms with an 
open d-shell. These values arc compared with the bind- 
ing energies calculated in Ref. 13[ and with semiempir- 
ical estimates of Ref. |14j. Within the accuracy of the 



current approach all results appear to be in good accord 
with each other. There is a correlation between the oc- 
cupation of the d-shell and the binding energies. In most 
cases the closer the occupation is to 100% the smaller is 
the binding energy, and the smaller arc the differences 
between the values obtained using different approaches. 
In the end, the results of the present work confirm the 
claim of Ref. [l3| that many open-shell atoms do bind 
the positron not only in the ground state but also in 
excited states. Note that according to our calculations, 
positrons do not bind to either Pd or Pt in their ground 
states. However, both atoms bind in excited states. 

The final energies of positron bound states and reso- 
nances are presented in Table Mil The resonant energy 
is the energy of the incident positron e for which capture 
into a bound state with an excited atom is possible. The 
resonant and binding energies are related to each other 
by 

e r = E CK - e b , (21) 

where E ex is the atomic excitation energy measured with 
respect to the ground state. To observe the resonance, 
one requires e r > 0, which means that the excitation en- 
ergy must be greater than binding energy (E cx > Sb)- In 
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TABLE II: Positron binding energies (BE, in eV) for open-shell atoms obtained using the SD equations (SD) with third-order 
correction (E3). In this table I is the ionization energy from the lowest state of a given configuration. The last collumn shows 
semiempirical values from Ref . [l4| . 



z 


Atom 


Valence 
configuration 


I 

(eV) 


„ a 
Oid 

(a.u.) 


SD 


BE, this work 
E3 


Total 


BE, other calculations 
Ref. [13] Ref. [14] 


21 


Sc 


3d4s 2 


6.56 


120 


0.908 


0.129 


1.037 




0.75(2) 






3d 2 4s 


5.13 




0.849 


0.109 


0.958 






22 


Ti 


3d 2 4s 2 


6.83 


98.5 


0.785 


0.110 


0.896 




0.84(3) 






3d 3 4s 


6.02 




0.727 


0.097 


0.825 






23 


V 


3d 3 4s 2 


6.74 


83.7 


0.678 


0.097 


0.775 




0.81(3) 






3d 4 4s 


6.48 




0.602 


0.086 


0.689 






24 


Cr 


3d 5 4s 


6.76 


78.3 


0.488 


0.077 


0.565 




0.54(8) 






3d 4 4s 2 


5.80 




0.590 


0.086 


0.676 






25 


Mn 


3d 5 4s 2 


7.43 


63 


0.496 


0.077 


0.574 




0.53(3) 






3d 6 4s 


5.32 




0.382 


0.068 


0.450 






26 


Fe 


3d 6 4s 2 


7.90 


56.7 


0.429 


0.069 


0.498 


0.28 


0.37(3) 






3d 7 4s 


7.04 




0.315 


0.061 


0.376 


0.09 




27 


Co 


3d 7 4s 2 


7.88 


50.7 


0.360 


0.061 


0.422 


0.26 


0.36(3) 






3d s 4s 


7.45 




0.243 


0.053 


0.297 


0.08 




28 


Ni 


3d 8 4s 2 


7.64 


45.9 


0.295 


0.055 


0.350 


0.24 


0.42(3) 






3d 9 4s 


7.61 




0.173 


0.046 


0.220 


0.07 




29 


Cu 


3d 10 4s 


7.72 


41 


0.125 


0.040 


0.166 


0.17(2) 


0.152(2) 


39 


Y 


4d5s 2 


6.22 


153 


0.845 


0.256 


1.102 




0.6(19) 






4d 2 5s 


4.86 




0.683 


0.258 


0.942 






40 


Zr 


4d 2 5s 2 


6.63 


121 


0.729 


0.209 


0.939 




0.8(11) 






4d 3 5s 


6.12 




0.623 


0.208 


0.831 






41 


Nb 


4d 4 5s 


6.76 


106 


0.527 


0.172 


0.699 




0.6(14) 






4d 3 5s 2 


6.62 




0.658 


0.178 


0.836 






42 


Mo 


4d 5 5s 


7.09 


86.4 


0.442 


0.145 


0.587 




0.45(4) 






4d 4 5s 2 


5.73 




0.583 


0.155 


0.739 






43 


Tc 


4d 5 5s 2 


7.28 


76.9 


0.461 


0.133 


0.594 


0.46 


0.62(4) 






4d 6 5s 


6.96 




0.355 


0.124 


0.479 


0.23 




44 


Ru 


4d 7 5s 


7.36 


65 


0.310 


0.109 


0.419 


0.21 


0.30(95) 






4d 6 5s 2 


6.43 




0.461 


0.121 


0.583 






45 


Rh 


4d 8 5s 


7.46 


58 


0.260 


0.095 


0.355 


0.20 


0.24(86) 


46 


Pel 


4d 10 


8.34 


32 


< 




< 




< 






4d 9 5s 


7.52 




0.205 


0.083 


0.288 










4d 8 5s 2 


5.23 




0.361 


0.097 


0.459 






71 


Lu 


5d6s 2 


5.43 


148 


0.188 


0.221 


0.410 




0.2(11) 


72 


Hf 


5d 2 6s 2 


6.83 


109 


0.305 


0.198 


0.503 




0.8(10) 






5d 3 6s 


5.08 




0.349 


0.190 


0.539 






73 


Ta 


5d 3 6s 2 


7.55 


88.4 


0.274 


0.166 


0.441 


0.45 


0.6(8) 






5d 4 6s 


6.34 




0.126 


0.147 


0.273 






74 


W 


5d 4 6s 2 


7.86 


74.9 


0.235 


0.141 


0.377 


0.46 


0.4(7) 






5d 5 6s 


7.50 




0.083 


0.118 


0.201 


0.30 




75 


Re 


5d 5 6s 2 


7.83 


65 


0.202 


0.121 


0.324 




0.42(12) 






5d 6 6s 


6.38 




0.051 


0.096 


0.147 






76 


Os 


5d 6 6s 2 


8.44 


57 


0.167 


0.105 


0.273 


0.47 


0.3(5) 






5d 7 6s 


7.80 




0.025 


0.079 


0.105 


0.29 




77 


Ir 


5d 7 6s 2 


8.97 


51 


0.137 


0.091 


0.229 


0.46 


0.2(4) 






5d 8 6s 


8.62 




0.005 


0.067 


0.072 


0.28 




78 


Pt 


5d 9 6s 


8.96 


44 


< 




< 


0.27 


< 






5d 8 6s 2 


8.86 




0.111 


0.080 


0.191 


0.46 








5d 10 


8.20 




0.020 


0.023 


0.044 


0.23 





"Ground-state atomic static dipole polarizabilities from Ref. 1 5 2f . 
''Configuration-interaction with many-body theory, Ref. [tJ. 
c Stochastic variation method, Ref. Q. 
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lighter atoms, the fine-structure splitting of the ground 
state is too small to satisfy this condition. However, in 
the heavier atoms, the fine structure can be larger than 
£6, and the corresponding low- lying excited states can 
form resonances. There is also another condition for ob- 
serving narrow resonances. The resonance energy should 
be smaller than the Ps formation threshold, s r < I— 6.8. 
For resonances lying above this energy, the Ps forma- 
tion channel will be open. This can reduce the resonance 
lifetimes and make the resonances too broad to be ob- 
served. Also, at these energies the positron-atom anni- 
hilation signal due to resonances will be "drowned" by 
the annihilation within the positronium. Table 11111 shows 
resonances which satisfy or approximately satisfy these 
two conditions. Note that we also show some states with 
negative e r values, which correspond to weakly bound 
positron states. Given the uncertainty in our calculation, 
they may in fact turn out to be low-lying resonances. 



open-shell atoms. The binding energies are in good agree- 
ment with previous estimates and indicate that atoms 
with open d-shells can bind positron not only in their 
ground but also in excited states. Many of the atoms 
considered appear to be good candidates for studying 
positron binding to atoms through resonant annihilation 
or resonant scattering. 
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TABLE III: Final binding and resonant energies (sb and e r , in eV) for the detection of positron-atom bound states through 
resonant annihilation or scattering. 





Atom 


Valence 


Excited states 


T7i a 

-c<ex 


£(. 




z 




configuration 










22 


Ti 


3(^4^ 


1 D 2 


0.90 


0.90 


0.0 






3d 3 4s 


5 F!, 5 F 2 


0.813, 0.818 


0.825 


-0.012, -0.007 


27 


Co 


3d 8 4s 


4p 4-p 4p 4r? 
^9/2, -T7/2j f 5/2i -T3/2 


0.43, 0.51, 0.58, 0.63 


0.30 


0.13, 0.21, 0.28, 0.33 








F7/2 


0.92 


0.30 


0.62 


28 


Ni 


3d 9 4s 




0.42 


0.22 


0.20 


43 


Tc 


4d 6 5s 


D5/2! D3/2, D1/2 


0.46, 0.50, 0.52 


0.48 


-0.02, 0.02, 0.04 


44 


Ru 


4d 7 5s 


5 F 1 


0.39 


0.42 


-0.03 








3 F 4 , 3 F 3 


0.81, 1.00 


0.42 


0.39, 0.58 








T 2 


1.00 


0.42 


0.58 


45 


Rh 


4d 8 5s 


F5/2, F3/2 


0.322, 0.43 


0.355 


-0.033, 0.08 








F7/2, F5/2 


0.71, 0.96 


0.36 


0.35, 0.60 


46 


Pd 


4d 9 5s 


2 [5/2] 3 , 2 [5/2] 2 


0.81, 0.96 


0.29 


0.52, 0.67 


72 


Hf 


5d 2 6s 2 


3 F 4 


0.57 


0.50 


0.07 


73 


Ta 


5d 3 6s 2 


F7/2, Fg/2 


0.49, 0.70 


0.44 


0.05, 0.26 








Pl/2, P3/2; P5/2 


0.75, 0.75, 1.15 


0.44 


0.31, 0.31, 0.71 








G7/2 

5 D 2 , 5 D 3 


1.20 


0.44 


0.76 


74 


W 


5d 4 6s 2 


0.41, 0.60 


0.38 


0.03, 0.22 








4 D 3 


0.77 


0.38 


0.39 








3 P2 


1.18 


0.38 


0.80 






5d 5 6s 


7 S 3 


0.36 


0.20 


0.16 


75 


Re 


5d 5 6s 2 


4 P 5 /2 


1.44 


0.32 


1.12 


76 


Os 


5d 6 6s 2 


5 D 2 , 5 D 3 , 5 D 1 , 5 Do 


0.34, 0.52, 0.71, 0.76 


0.27 


0.07, 0.25, 0.44, 0.49 








3 7 t 3tt 3tt 


1.78, 1.84, 1.84 


0.27 


1.51, 1.57, 1.57 






5d 7 6s 


5 F 5 , 5 F 4 


0.64, 1.08 


0.10 


0.54, 0.98 


77 


Ir 


5d 7 6s 2 


F3/2, F5/2, F7/2 


0.51, 0.72, 0.78 


0.23 


0.28, 0.49, 0.55 








2 /~i 2 

^9/2] W/2 


1.73, 2.20 


0.23 


1.50, 1.97 








4p 4p 4p 

r 5/2i r 3/2j -Tl/2 


2.00, 2.30, 2.51 


0.23 


1.77, 2.07, 2.27 








Hn/2 


2.43 


0.23 


2.20 






5d 8 6s 


4p 4-p 4p 4p 
-T9/2! -T7/2; -T5/2, ^3/2 


0.35, 0.88, 1.22, 1.47 


0.07 


0.28, 0.81, 1.15, 1.40 








2p 2p 
r 3/2i r l/2 


1.31, 1.55 


0.07 


1.24, 1.48 








2 P 2tti 

-T5/2, -T7/2 


1.51, 1.62 


0.07 


1.44, 1.55 








P5/2 


1.60 


0.07 


1.53 


78 


Ft 


5d 8 6s 2 


3p 3tt> 


1.25, 1.92 


0.19 


1.06, 1.73 








3 P 2 


0.81 


0.19 


0.62 






5d 10 




0.76 


0.04 


0.72 



"Atomic excitation energy with respect to the ground state from 
Ref. [H. 
^Resonance energy from Eq. 121> . 
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